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Summary 

t  he  authors  report  an  interdecadal  oscillation  in  a  wind-  and  thermally-driven  ocean  general  circulation  model 
(OGCM).  The  oscillation  ix  tantalizing  in  that  it  occurs  under  a  relatively  strong  thermal  damping  (26.3  WmA"1) 
Examinations  involving  a  two-dimensional  OGCM.  a  simple  thermal  ' flip-flop '  model,  and  a  three-dimensional 
OGCM  with  and  without  the  nonlinear  effect  of  temperature  in  the  state  equation  of  sea  water  demonstrate  that  the 
oscillation  is  not  driven  h\  mechanisms  such  as  the  so-called  convective  oscillator,  or  the  advectivc  overshooting 
oscillator.  Instead,  the  oscillation  is  associated  with  the  propagation  of  modelled  long  RosSby-waves  with  westward 
phase  velocity.  It  is  fount!  that  the  north  -south  basin  extent  of  the  OGCM  is  an  important  factor  a  larger  north- 
south  basin  extent  is  conducive  to  the  generation  of  thermal  wind  flows  crucial  for  initiating  the  propagation  of 
disturbances.  The  geophysical  relevance  of  these  results  is  discussed. 
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1 .  Introduction 

Interdecadal  variabilities  have  been  studied  extensively  in  ocean  general  circulation 
models  (OGCMs).  Oscillations  induced  by  saline  forcing  (e.g.  Marotzke  1990;  Weaver 
and  Sarachik  1991a.b:  Wright  and  Stocker  *1991;  Winton  and  Sarachik  1993;  Huang  1994; 
Huang  and  Chou  1994;  Cai  1995;  Chen  and  Ghil  1995)  are  usually  examined  in  models 
under  so-called  mixed  boundary-conditions  (that  is.  the  use  of  a  restoring  boundarv- 
condition  on  surface  temperature  and  a  flux  boundary -condition  on  surface  salinity  (Bryan 
1986)).  Oscillations  induced  by  thermal  forcing  (e.g.  Cai  et  al.  1995;  Greatbatch  and  Zhang 
1995:  Cai  1996)  have  been  examined  in  a  purely  thermal,  or  thermally-dominant,  system. 
Most  of  these  studies  explain  the  oscillation  as  driven  by  an  advectivc  overshooting  oscil¬ 
lator.  In  a  thermal  system,  the  argument  goes  as  follows.  When  the  overturning  is  stronger 
than  the  average,  heat  storage  increases  al  high  latitudes.  This  reduces  the  contrast  be¬ 
tween  densities  at  pole  and  equator,  and  weakens  the  overturning  circulation,  which,  in 
turn,  reduces  the  poleward  heat-transport  and  leads  to  a  high-latitude  cooling.  The  pole 
to  equator  density  contrast  then  increases,  accelerating  the  overturning,  and  the  process 
repeats  itself. 

Lenderink  and  Haarsma  (1994)  and  Pierce  et  al.  (1995)  considered  the  competing 
roles  of  thermal  and  saline  forcing  components  in  determining  the  stability  of  a  water  col¬ 
umn.  Both  studies  emphasized  the  role  of  a  convective  oscillator.  Lenderink  and  Haarsma 
(1994)  found  that,  depending  upon  the  comparative  role  of  the  two  forcing  components,  the 
water  column  could  have  three  different  regimes:  one  permanently  convective,  one  perma¬ 
nently  non-convective,  and  one  oscillatory  (periodically  switching  between  the  convective 
and  non-convective  regimes).  Pierce  et  al.  (1995)  proposed  another  kind  of  convective 
oscillator.  They  showed  that,  in  order  for  oscillations  to  occur  in  their  model,  the  nonlin¬ 
ear  effect  of  temperature  in  the  state  equation  of  sea  water  is  a  necessary  condition.  The 
mechanism  is  supported  by  an  analytical,  two-layered  ‘flip-flop'  model. 

Recently,  the  role  of  the  coastal  wave-guide  has  been  emphasized.  This  is  associated 
w  ith  the  thermohalinc  adjustment  in  the  generation  of  oceanic  variability  (Greatbatch  and 
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Peterson  1996:  Winton  1996).  The  adjustment  problem  was  first  studied  by  Davey  (1983) 
and  further  explored  by  Wajscwicz  and  Gill  (1986).  In  a  spin-down  of  an  OGCM  driven  by 
a  density  field  that  varies  only  in  the  meridional  direction,  Wajsowicz  and  Gill  (1986)  found 
the  importance  of  viscous,  baroclinic  boundary-waves  in  the  model  adjustment  process. 
The  problem  has  also  been  addressed  by  Winton  (1996).  He  extended  the  early  studies 
to  include  the  the  propagation  of  boundary  waves  in  regions  of  formation  of  weakly-  or 
non-stratified  deep-water  in  models  driven  by  a  constant  heat-flux  forcing  alone.  He  found 
that  under  a  constant  heat-flux  forcing,  generation  of  interdecal  oscillation  was  a  robust 
feature.  It  involves  periodic  growth  and  decay  of  a  baroclinic  disturbance,  which  propagates 
cyclonically  around  the  basin.  The  result  is  confirmed  by  Grcatbatch  and  Peterson  (1996) 
in  an  OGCM.  again  under  a  constant  heat-flux  forcing.  These  authors  believed  that  the 
waves  in  their  studies  were  Kelvin  waves,  although  strictly  speaking  Kelvin  waves  were 
not  resolved  in  their  models  because  of  the  low  resolution. 

Griffies  and  Tziperman  (1995)  interpreted  the  oscillation  in  the  Geophysical  Fluid 
Dynamics  Laboratory  (GFDL)  coupled  atmosphere-ocean  model  (Delworth  el  al.  1 993)  as 
an  external  excitation  of  a  damped  oscillatory  thermohaline  mode  by  the  stochastic  effects 
of  the  atmospheric  variability.  This  is  in  sharp  contrast  to  the  explanation  given  by  Cai 
el  al.  (1995)  and  Greatbatch  and  Zhang  (1995).  who  concluded  that  the  oscillation  might 
be  an  internal  oceanic  process.  In  their  ocean-onlv  models,  driven  by  a  constant  heat-flux 
forcing,  they  were  able  to  produce  many  oscillation  features  similar  to  those  observed  in 
the  North  Atlantic  (Deser  and  Blackmon  1993;  Kushnir  1994).  and  to  those  of  the  inter- 
decadal  oscillations  in  the  GFDL  coupled  model  (Delworth  et  al.  1993).  Latif  and  Barnett 
(1994.  1996)  focused  on  interdecadal  climate  variability  over  the  North  Pacific  Ocean  and 
North  America.  They  found  that  the  modes  involve  unstable  ocean-atmosphere  interactions 
over  the  North  Pacific,  and  that  the  time  scale  of  the  oscillation  is  mainly  determined 
by  that  of  baroclinic  planetary  waves,  which  are  generated  by  anomalies  in  w  ind-stress 
curl. 

In  this  study,  we  report  an  interdecadal  oscillation  under  a  restorative  forcing  in  an 
OGCM  similar  to  that  reported  by  Cai  et  al.  (1995).  The  tantalizing  aspect  of  the  oscillation 
is  that  it  occurs  under  a  relatively  strong  thermal  damping  (26.3  W  m':K_1).  This  result 
is  in  sharp  contrast  to  those  found  in  previous  studies  (e.g.  Cai  et  al.  1995;  Greatbatch  and 
Peterson  1996:  Winton  1996).  In  these  studies,  a  flux  condition  was  required  for  persistent 
oscillations  and  the  possibility  of  oscillations  under  a  restorative  forcing  was  not  examined. 

As  can  be  seen  from  the  above  discussion,  there  are  a  number  of  possible  mecha¬ 
nisms  for  oscillations.  Which  of  the  above,  if  any.  is  responsible  for  the  generation  of  the 
oscillation  in  the  present  study?  An  important  clue  is  that  oscillations  arc  easily  induced 
in  the  present  model  but  were  not  so  in  the  work  reported  by  Cai  et  al.  (1995).  We  begin 
our  exploration  by  describing  the  difference  in  model  configurations  between  the  present 
study  and  that  of  Cai  et  al.  (1995)  in  section  2.  In  section  3.  we  show'  the  generation  of  the 
oscillation  under  a  relatively  strong  thermal  damping  rate,  and  the  low  dependence  upon 
initial  conditions.  In  section  4,  possible  causes  are  examined.  We  show  that  the  oscillation 
in  our  study  is  not  driven  by  a  convective  oscillator,  nor  by  an  advectivc  overshooting  os¬ 
cillator.  In  section  5.  we  demonstrate  that  the  oscillation  is  associated  with  boundary-wave 
propagation,  and  that  the  wave  is  a  long  Rossby-wave  (with  w'estward  phase  velocity), 
rather  than  a  Kelvin  wave  as  in  the  studies  of  Winton  (1996)  and  Greatbatch  and  Peterson 
(1996).  We  show'  that  basin  length  is  an  important  factor  for  the  generation  of  the  oscil¬ 
lation.  A  larger  north-south  extent  of  the  ocean  basin  enables  this  mechanism  to  operate 
so  effectively  that  it  can  sustain  a  considerable  damping,  and  that  a  flux  condition  is  not 
required.  The  implication  and  the  geophysical  relevances  of  the  oscillation  are  discussed 
in  section  6. 
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TABLE  1.  Distribution  of  vertical 
levels 


Level 

Thickness 

(m) 

Depth  of  (T,  S ) 
(m) 

1 

25.0 

12.5 

2 

25.0 

37.5 

3 

40.0 

70.0 

4 

70.0 

125.0 

5 

1 10.0 

215.0 

6 

200.0 

370.0 

7 

330.0 

635.0 

8 

450.0 

1025.0 

9 

650.0 

1575.0 

10 

900.0 

2350.0 

11 

900.0 

3250.0 

12 

900.0 

4150.0 

TABLE  2.  Standard  values  of  model  COEF¬ 
FICIENTS 


Value 

Parameter 

Symbol 

(nr  s'1) 

Horizontal  diffusivity 

Ath 

1  x  10-1 

Horizontal  viscosity 

Amh 

3  x  10s 

Vertical  diffusivity 

Atv 

1  x  10"* 

Vertical  viscosity 

-4  mV 

I  x  lO-1 

2.  The  model 

This  study  employs  the  version  of  the  Bryan-Cox-Semtner  OGCM  used  by 
Pacanowski  el  at  (1991 ).  which  is  based  on  the  work  of  Bryan  (1969).  The  state  equation 
is  that  described  by  Gill  (1982),  and  includes  the  full  nonlinear  effect  of  temperature.  In  all 
experiments,  there  is  no  salinity  forcing  and  the  salinity  is  kept  constant  and  uniform.  The 
horizontal  grid  spacing  is  4  of  latitude  by  4  of  longitude.  The  model  has  twelve  levels  in 
the  vertical,  at  depths  listed  in  Table  1.  Values  assigned  to  the  various  model  parameters  are 
listed  in  Table  2.  These  are  the  same  as  those  used  by  Cai  (1995),  Cai  and  Godfrey  (1995) 
and  Cai  eta!.  (1995).  The  model  domain  is  a  flat-bottomed,  two-hemisphere  Atlantic  basin. 
It  has  a  width  of  60  and,  in  the  basic  cases,  a  length  of  144  extending  from  72CS  to  72UN. 
This  extends  the  model  used  by  Cai  etui.  (1995)  in  the  north-south  direction  by  one  model 
grid-point  in  each  hemisphere.  In  some  cases  the  basin  length  is  shortened  or  extended 
for  testing  the  sensitivity  to  basin  length.  The  southern  hemisphere  includes  a  modelled 
Antarctic  Circumpolar  Current  (ACC)  passage  from  48  S  to  64  S.  and  a  sill  2350  m  deep 
in  the  model's  Drake  Passage.  Cyclic  conditions  are  applied  there.  This  sill  is  intended  to 
provide  a  realistic  ACC  transport  by  setting  up  a  bottom  pressure  difference  between  the 
two  sides  of  the  passage  south  of  Cape  Horn  (Gill  and  Bryan  1971;  Holland  1973;  Cai 
1994).  Otherwise,  the  flat  bottom  would  lead  to  an  unrealistically  large  wind-driven  ACC: 
so  one  has  to  specify  an  ACC  (see  Cai  and  Godfrey  1995).  Another  small  modification 
from  the  model  used  by  Cai  et  al.  (1995)  is  the  inclusion  of  a  model  Weddell  Sea  to  allow 
an  Antarctic  Circumpolar  Convective  cell.  (Detailed  layout  can  be  seen  in  Fig.  7.) 
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Figure  I  (a)  Surface  relaxation  temperature  (  C).  and  (b)  overturning  stream  function  (Sv)  for  R!5o.  that  is. 
under  a  restoration  with  a  restoring  time  of  15  days,  and  initiated  from  a  resting  ocean. 


TABLE  3.  Details  of  model  experiments  discussed  in  section  3 


Run 

Initial 

Condition 

Model 

Domain 

Forcing 

Persistent 

Oscillations 

RI5o 

Rest 

72S  to  72 'N 

Restoring  (15  d) 

none 

R45 

R15o 

72  S  to  72  N 

Restoring  (45  d) 

24.7  \ 

R45i> 

Rest 

72  S  to  72  N 

Restoring  (45  d) 

23.7  v 

R30o 

Rest 

72  S  to  72  N 

Restoring  (3(1  d) 

none 

R45oShtB 

Rest 

64  S  to  64  N 

Restoring  (45  d) 

none 

2DR45o 

Rest 

72  S  to  72  N 

Restorine  (45  d) 

none 

LnrR43ci 

Rest 

72  S  to  72  N 

Restoring  (45  d) 

24.6  y 

R45oflt 

Rest 

72  S  to  72  N 

Restoring  (45  d) 

none 

R45ShtBp 

R45oShlB 

64  S  to  64  N 

Restoring  (45  d) 

none 

R3()p 

R30o 

72S  to  72  N 

Restoring  (30  d) 

none 

R45o0.5u> 

Rest 

72  S  to  72  N 

Restoring  (45  d) 

16.8  v 

R45ol.5u> 

Rest 

72  S  to  72  N 

Restoring  (45  d) 

37  Jy 

The  fifth  column  indicates  whether  or  not  persistent  oscillations  occurred.  If  yes. 
the  period  is  given  in  years  (y). 


3.  Initial  model  results 

(a)  Model  spin-up 

Initially,  the  OGCM  was  spun  up  by  restoring  the  temperature  at  the  uppermost  level 
to  a  zonaily  Uniterm  value  given  by  the  profile  shown  in  Fig.  1(a).  which  approximates 
the  observed  zonaily  averaged  field.  This  run  is  referred  to  as  R15o  (see  Table  3),  and 
starts  from  a  resting  ocean  (denoted  by  ‘o  ).  The  profile  has  the  analytical  expression 
used  by  Cai  et  al.  (1995)  and  the  expression  applies  to  the  extended  model  domain  from 
68  to  72\  The  restoring  time  is  15  days  (hence  the  '15'  in  R15o).  This  corresponds 
to  a  thermal  coupling  strength  of  79  W  nr2K_1.  The  model  is  also  subject  to  the  wind 
stress  used  by  Bryan  (1987).  The  wand  stress  is  symmetric  about  the  equator,  and  mimics 
well  the  observed  wind  stress  with  features  of  strong  westerlies  at  mid-to-high  latitudes 
and  moderate  easterlies  in  the  tropics.  In  all  model  runs,  the  wind  is  kept  constant.  The 
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lower-level  acceleration  techniques  of  Bryan  ( 1984)  are  not  used  at  any  stage.  After  about 
4000  years  of  integration,  the  spin-up  reaches  a  statistically  steady  state.  The  overturning 
circulation  is  shown  in  Fig.  1(b),  and  is  qualitatively  similar  to  that  in  the  real  Atlantic 
Ocean.  It  features  a  northern  sinking  cell,  an  ACC  convective  cell,  and  a  northern  intrusion 
cell  from  the  south.  The  model  ACC  reaches  104  Sv.  The  modelled  circulation  is  of  the 
conveyor-belt  type  and  is  similar  to  that  described  by  Cai  et  al.  (1995). 

(h)  Oscillation  under  u  restorative  condition 

Initiated  from  the  steady  state  of  R I  So,  R45  (Table  3)  was  carried  out.  in  which  the 
temperature  at  the  uppermost  level  was  relaxed  to  the  same  climatology  (Fig.  1(a))  but  with 
a  restoring  time  of  45  days.  At  the  moment  of  changing  the  surface  forcing,  the  surface 
heat-flux  was  suddenly  and  uniformly  reduced  to  one  third  of  that  before  the  change. 
The  reduction  of  the  heat  flux  shocked  the  system  somewhat  initially  but.  as  it  recovered, 
oscillations  with  a  period  of  24.7  years  appeared.  These  can  be  seen  in  Fig.  2(a).  which 
shows  the  time  scries  of  the  overturning  sampled  at  a  location  (60  N.  2350  m)  near  the 
place  where  the  overturning  of  the  steady  state  of  R15o  is  a  maximum.  Unless  otherwise 
stated,  all  single-point  time-series  of  overturning  were  taken  at  this  location. 

Since  R45  was  initiated  from  the  steady  state  of  R15o.  a  question  arises  as  to  what 
would  happen  if  the  OGCM  were  spun  up  from  the  same  steady  state  as  R15o.  that  is  a 
resting  ocean.  To  this  end.  R45o  (Table  3)  was  carried  out.  Figure  2(b)  shows  the  time  series 
of  the  overturning.  Oscillations  persist.  The  period  is  23.7  years.  Thus  the  oscillation  is 
almost  independent  of  the  initial  condition.  In  fact,  from  Fig.  2(b)  we  sec  that  the  oscillation 
occurred  at  the  early  stage  of  the  spin-up  and  initially  had  large  amplitudes.  The  amplitudes 
reduced  as  the  integration  proceeded  until  a  stage  was  reached  when  they  stablized.  The 
reduction  was  due  to  the  weakening  (damping)  effect  associated  with  the  restoring  forcing. 
This  can  be  further  illustrated  from  a  time  series  (Fig.  3(a))  for  run  R30o  (Table  3).  This  run 
was  identical  to  R45o.  except  that  the  restoring  time  was  30  days.  Initially  the  evolution 
of  the  overturning  resembled  that  of  R45o  but  the  oscillation  was  eventually  damped. 


Tim  (yeAra I  7c»*>  (years) 

Figure  2.  (a)  Time  series  of  the  overturning  (Sv)  for  run  R45.  that  is.  initiated  from  the  steady  state  of  RI5o  and 
under  a  restoration  with  a  restoring  time  of  45  days.  The  overturning  is  sampled  at  a  location  (60  N.  2.730  m)  near 
the  place  where  the  overturning  of  the  steady  state  of  R 1 5o  is  a  maximum,  (b)  Time  scries  of  the  overturning  (Sv) 
for  run  R45o,  that  is.  under  a  restoration  with  a  restoring  time  of  45  days  but  initiated  from  a  resting  ocean. 
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Figure  3.  (a)  Time  series  of  ihe  overturning  (S'  )  for  R30o.  that  is  under  a  restoration  wilh  a  restoring  lime  of  30 

days  and  initiated  from  a  resting  ocean,  (b)  Time  series  of  the  overturning  for  run  R45oShlB.  that  is  in  a  model 
with  a  short  north-south  basin  extent,  under  a  restoration  with  a  restoring  time  of  45  days,  and  initiated  from  a 

resting  ocean. 


The  oscillation  is  interesting,  especially  when  compared  with  the  result  of  Cai  et  al. 
(1995)  who  found  that,  even  under  a  diagnosed  heat-flux  forcing  (without  the  damping 
effect  from  the  surface  restoration)  no  oscillation  occurred.  A  heat-flux  redistribution  w  as 
necessary  for  the  generation  of  oscillations.  As  commented  upon  earlier,  the  only  difference 
between  the  present  model  and  that  of  Cai  et  al.  (1995  )  is  that  the  north-south  extent  of 
the  model  basin  is  4  (one  grid-length)  longer  in  the  present  model  in  each  hemisphere. 
The  extent  of  model  basin  appears  to  be  important.  To  test  this,  run  R45oShtB  (Tabic  3) 
was  carried  out.  This  run  was  identical  to  R45o  except  that  the  basin  extent  was  from  64  S 
to  64  N,  two  grid-lengths  shorter  than  that  in  the  R45o  case.  This  is  chosen  to  facilitate 
a  clear-cut  comparison  between  the  short  basin  (R45oShtB)  and  the  long  basin  (R45o) 
cases.  Since  the  restoring  temperature  took  the  same  profile  as  in  R45o  (Fig.  1(a)).  the 
restoring  temperature  at  the  polar  boundary  in  this  run  was  warmer  than  that  in  R45o. 
Figure  3(b)  shows  the  evolution  of  the  overturning.  We  sec  that  no  persistent  oscillation 
was  produced.  Note  that  the  maximum  overturning  in  the  northern  hemisphere  in  R45oShtB 
is  of  a  comparable  size  to  that  in  R45o,  but  moves  southward,  so  that  the  overturning  shown 
is  smaller. 


4.  Some  possible  causes  of  the  oscillation 

Consider  run  R45o  and  R45oShtB.  One  may  ask  several  questions.  Is  it  the  lower 
restoring  temperature  at  the  polar  boundary,  or  the  associated  change  in  wind,  that  is 
important  as  a  result  of  the  north-south  extension  of  the  basin?  Or  is  it  that  the  exact 
latitude  is  important?  Or  is  it  rather  that  a  large  initial  extent  of  the  convective  region 
predisposes  the  model  towards  increased  variability,  since  convection  is  intermittent  even 
at  the  spin-up  equilibrium?  In  order  to  answer  these  questions,  many  experiments  have  been 
carried  out  using  both  the  long-  and  short-basin  configurations  of  the  model.  The  results 
indicate  that  it  is  the  extent  of  the  basin,  rather  than  changes  in  the  restoring  temperature 
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or  the  wind  stress,  which  accounLs  for  the  different  behaviour  of  the  circulation  between 
the  short-  and  long-basin  cases. 

A  significant  contrast  in  the  evolution  of  the  overturning  between  long-basin  cases 
(e.g.  R45o)  and  short-basin  cases  (e.g.  R45oShtB')  is  that  the  oscillations  develop  in  the 
initial  stage  (e.g.  the  first  500  years)  in  the  former  but  are  absent  from  the  early  stage  of 
integration  in  the  latter.  This  suggests  that,  from  the  initial  stage,  the  adjustment  process 
in  the  long-  and  short-basin  cases  is  quite  different. 

In  what  follows,  we  show  that  the  oscillation  in  our  study  is  not  driven  by  a  convective 
oscillator,  nor  by  an  advcctivc  overshooting  oscillator.  Since  we  are  dealing  with  a  thermal 
system,  the  convective  oscillator  proposed  by  Lenderink  and  Haarsma  (1994)  does  not 
operate  here.  Work  by  Pierce  et  ul.  (1995)  using  a  'flip-flop’  mode!  implicitly  suggested 
that,  when  density  is  a  nonlinear  function  of  temperature,  oscillations  may  be  induced,  even 
in  a  purely  thermal  circulation.  1  n  all  the  experiments  described  so  far,  the  nonlinear  effect  of 
temperature  in  the  state  equation  is  included.  In  order  to  test  if  this  oscillator  is  responsible 
for  the  oscillation  in  our  OGCM,  an  analysis  was  carried  out  using  a  simplified  thermal- 
only  version  of  the  two-layered  flip-flop  model.  Details  are  given  in  an  appendix.  The 
simplified  flip-flop  model  allowed  a  surface  cooling  and  the  associated  heat  loss.  In  order 
to  maintain  a  heat  budget,  the  lost  heat  was  re-introduced  into  the  system  as  a  subsurface 
healing.  As  shown  in  the  appendix,  because  of  the  nonlinear  effect  of  the  temperature  in  the 
state  equation,  such  a  thermal  system  may,  depending  on  the  magnitude  of  the  subsurface 
healing,  have  three  different  regimes:  one  permanently  convective,  one  permanently  non- 
convective,  and  one  oscillatory  (periodically  switching  between  the  convective  and  non- 
convective  regimes). 

There  are  at  least  two  ways  to  check  if  the  mechanism  works  in  our  OGCM.  Firstly, 
the  flip-flop  model  suggests  that  once  the  nonlinear  effect  of  temperature  in  the  state 
equation  is  removed,  no  oscillation  occurs.  However,  when  R45o  was  repeated  with  a  den¬ 
sity  equation  set  to  that  in  the  flip-flop  model,  giving  run  LnrR45o  (Table  3).  oscillations 
still  occurred.  This  indicates  that  the  nonlinearity  is  not  what  drives  the  oscillation  in  our 
OGCM.  Additional  support  can  be  found  from  the  work  of  Huang  and  Chou  (1994)  who 
used  a  salinity-driven  model,  in  which  persistent  oscillations  were  produced  despite  the 
state  equation’s  being  nearly  linear  in  salinity.  Secondly,  if  the  mechanism  in  the  flip-flop 
model  is  what  drives  the  oscillations  in  our  OGCM.  then  oscillations  should  also  occur  in  a 
two-dimensional  version  of  the  OGCM.  To  test  this  inference,  we  carried  out  run  2dR45o 
(Table  3).  The  cast-west  extent  was  reduced  to  two  model  tracer-grids  (that  is,  only  one 
velocity  grid)  and  there  were  no  east  or  west  boundaries.  Other  details  were  identical  to 
R45o.  No  persistent  oscillation  resulted,  further  confirming  that  the  convective  oscillator 
displayed  in  the  flip-flop  model  was  not  the  cause  of  the  oscillation  in  our  OGCM.  Indeed, 
an  examination  shows  that  the  subsurface  heating  obtained  from  the  OGCM  near  the  polar 
regions  was  generally  greater  than  that  required  for  oscillations  to  occur  in  the  flip-flop 
model.  The  no-oscillation  solution  in  2dR45o  also  suggests  that  the  advective-overshooting 
oscillator  is  not  what  drives  the  oscillation  in  our  OGCM  cither. 

5.  Oscillation  and  wave  propagation 

A  detailed  examination  reveals  that  the  oscillation  in  the  present  study  is  associated 
with  waves  that  propagate  along  the  model  northern  boundary.  The  wave  propagation  for 
run  R45  can  be  seen  in  Fig.  4.  It  shows  the  evolution  of  the  rates  of  change  of  heat  content 
HC,  which  is  defined  as 
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Figure  4.  £  volution  of  the  rates  of  change  of  heat  content  over  water  columns  along  water  grids  nearest  to  the 

western,  northern  and  eastern  boundaries.  The  distance  is  measured  clockwise  from  the  south-western  end  of  the 
model  domain  and  is  plotted  as  3  straight  line.  The  contour  interval  is  S  W' m-2. 


Here  cp  is  the  specific  heal  capacity  and  At  is  the  reference  density.  Rates  of  change  in 
heal  content  are  therefore  described  by  9HC/3r.  Shown  in  Fig.  4.  is  this  quantity  over 
water  columns  at  grid  points  nearest  to  the  western,  northern,  and  eastern  land  boundaries. 
The  distance  along  the  boundaries  is  measured  clockw  ise  from  the  south-western  end  of 
the  model  domain  and  is  plotted  as  a  straight  line.  The  plot  covers  a  period  from  vear 
3545  to  year  3577  (see  Fig.  2(a)).  Several  features  emerge.  Firstly,  most  of  the  time  of 
an  oscillation  is  spent  in  propagation  along  the  weakly  or  non-stratified  boundary,  that 
is  along  the  northern  boundary.  Along  this  boundary,  the  propagation  speed  is  slow'  and 
the  amplitude  large.  By  contrast,  along  stratified  boundaries  the  propagation  speed  is 
much  faster  and  the  amplitude  much  smaller.  Secondly,  anomalies  propagate  faster  during 
increasing  heat-content  than  during  decreasing  heat-content.  Finally,  an  oscillation  period 
is  the  sum  of  the  time  for  positive  and  negative  anomalies  to  cross  the  basin. 

To  reveal  the  interrelationship  between  the  wave  properties  and  the  thcrmohaline 
circulation,  the  evolution  of  the  overturning  covering  the  same  period  is  shown  in  Fig. 
5(a).  and  the  anomalies  of  overturning  and  sea-surface  temperature  (SST)  from  a  mean 
circulation  every  quarter  period  arc  shown  in  Figs.  6  and  7  respectively.  The  overturning 
reaches  a  maximum  at  year  3550  and  a  minimum  at  year  3560.  A  quarter  period  before  the 
overturning  reaches  the  maximum,  the  northern  boundary'  is  occupied  by  cool  anomalies 
(as  can  be  inferred  from  Fig.  7(d)).  This  means  that,  at  that  time,  the  anomaly  of  the  north- 
soulh  pressure  gradient  is  largest.  From  then  until  the  time  of  maximum  overturning,  the 
meridional  northward  flow  increases  (Fig.  6(d))  and  this  is  associated  with  an  increase  in 
the  east-west  pressure  gradient  and  the  development  of  a  positive  SST  anomaly  in  the  east 
(Fig.  7(d))  on  the  way.  The  anomalous  northward  flow  not  only  pumps  a  heat  anomaly  into 
the  north-eastern  corner,  but  also  causes  the  anomaly  to  propagate  westwards  (Fig.  4).  The 
weak  or  non-existent  stratification  along  the  northern  boundary  impedes  the  movement 
along  the  leading  edge,  allow  ing  more  heal  to  be  accumulated  at  the  front,  leading  to  the 
increase  in  amplitude.  Then,  from  year  3550  to  year  3554,  the  heat  content  increases.  At 
year  3552  and  near  the  north-western  aimer,  it  increases  at  the  greatest  rate.  An  animation 
clearly  shows  that  the  wave  front  turns  southward  in  this  corner.  It  then  propagates  along 
the  western  boundary  with  increasing  speed  and  reducing  amplitude.  During  this  time, 
the  northern  boundary  is  occupied  by  positive  SST  anomalies  (Fig.  7(a.  b)).  The  large 
heat  increase  in  the  north-western  corner  then  leads  to  a  decrease  in  the  east-west  pressure 
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Figure  5.  (a)  Time  series  of  Ihe  overturning  for  R45  covering  a  33-year  period  starling  from  year  3545.  The 
overturning  reaches  a  maximum  ar  year  3550  and  a  minimum  at  year  3560.  (h)  Time  series  of  the  temperature 
averaged  over  the  northern  half  of  the  northern-hemisphere  model  basin,  subtracted  from  the  temperature  over  the 
southern  half  of  the  northern-hemisphere  model  basin.  The  maximum  of  the  associated  density -gradient  leads  the 
maximum  of  the  overturning  by  about  6  years,  that  is  by  about  90  . 


gradient  and  in  the  overturning  (Fig.  6(a,  b)).  From  year  3554  to  year  3560.  the  overturning 
decreases  from  below  the  mean  (Fig.  6(b,  d))  to  a  minimum,  and  cold  anomalies  develop 
at  the  north-eastern  corner  and  propagate  westwards,  associated  with  the  opposite  phase 
of  the  oscillation. 

The  feature  of  wave  propagation  can  also  be  seen  in  the  phase  relationship  between 
the  strength  of  the  overturning  circulation  and  the  north-south  density  difference.  Figure 
5(b)  shows  the  difference  between  temperatures  averaged  over  the  northern  and  southern 
halves  of  the  northern-hemisphere  model  basin.  Comparing  panels  (a)  and  (b)  of  Fig.  5, 
we  see  that  the  maximum  in  the  density  gradient  leads  the  maximum  overturning  by  about 
6  years,  that  is  by  about  90  . 

The  above  features  are  somewhat  similar  to  those  found  bv  Winton  (1996)  and  Great- 
batch  and  Peterson  (1996).  However,  the  wave  in  these  two  studies  was  a  Kelvin  wave.  In 
order  to  identify  the  nature  of  the  wave  in  the  present  study.  R45of0  was  carried  out.  This 
experiment  was  identical  to  R45o  except  that,  north  of  22  N,  /{-effects  were  switched  off, 
and  a  constant  Coriolis  parameter  (that  for  22  N)  was  used.  No  oscillation  was  generated. 
The  generation  of  Kelvin  waves  does  not  require  the  presence  of  /i -effects.  The  lack  of 
wave  generation  when  /{-effects  were  absent  suggests  that  the  wave  in  the  present  study 
was  not  a  Kelvin  wave.  Instead  it  was  a  long  Rossby  wave  with  westward  phase  velocity. 
In  the  absence  of  any  bottom  topography,  such  long  Rossby  waves,  may  be  generated  in 
the  presence  of  /{-effects.  This  is  so  in  our  OGCM. 

The  Rossby  wave  in  the  present  study  was  forced  by  density  anomalies  in  the  north¬ 
eastern  corner  of  the  basin.  This  being  so.  the  speed  of  anomalies  across  the  basin  was 
determined  by  the  slowly  varying  forcing  and  did  not  need  to  match  that  of  a  free  Rossby 
wave.  However,  in  order  to  gain  an  impression  of  the  relationship  between  the  oscillation 
period  and  the  phase  velocity  of  the  model  Rossby  waves,  we  approximate  the  phase 
velocity  of  the  modelled  waves  by  that  of  a  free  Rossby  wave.  For  a  free  Rossby  wave, 
the  phase  velocity  may  be  expressed  in  terms  of  /Jo2/ fa \  where  c  is  the  stratified  gravity 
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Figure  6.  Anomalies  of  overturning  (Sv  )  from  a  mean  circulation  for  R45  at  about  a  quarter  period  interval:  (a) 
at  year  3550:  (b)  at  year  3554;  (c)  at  year  3560;  (d)  at  year  3560.  Panels  (a)  and  (d)  show  the  anomalies  when  the 
overturning  is  at  a  maximum  and  a  minimum  respectively. 


wave  speed  and  /0  is  the  Coriolis  parameter  (Gili  1982).  Since  anomalies  propagate  faster 
during  increasing  heat-content  than  during  decreasing  heat-content,  and  since  an  oscilla¬ 
tion  period  is  the  sum  of  the  time  for  positive  and  negative  anomalies  to  cross  the  basin,  we 
obtained  from  R15  two  vertical  profiles  of  zonal-mean  density  of  the  northern-most  grids, 
one  averaged  over  the  period  of  increasing  heat  content  and  one  over  the  period  of 
decreasing  heat-content.  The  two  profiles  were  then  used  to  calculate  the  stratified  gravity- 
wave  speed  during  increasing  (n)  and  decreasing  (cD)  heat-content.  We  then  employed 
the  method  for  decomposition  of  vertical  modes  detailed  by  Gill  (1982)  and  boundary 
conditions  applicable  to  the  rigid  lid  and  a  flat  ocean-bottom  (Moore  and  Philander  1977). 
The  two  profiles  give  the  first  mode  C|  and  Cp  of  0.441  and  0.392  ms-1,  corresponding  to 
times  needed  to  travel  from  east  to  west  along  the  northern  boundary,  using  a  latitude  0  of 
70' ,  of  9.03  and  11.44  years,  respectively.  The  total  lime  is  therefore  20.7  years,  close  to 
the  oscillation  period. 
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The  oscillation  mechanism  is  clearly  quite  sensitive  to  the  presence  of  fi  effects.  In 
order  to  ensure  that  the  zero-beta  change  in  run  R450R)  did  not  have  influences  beyond 
suppressing  the  Rossby  waves  and  to  further  confirm  the  Rossby-wave  mechanism,  we 
have  run  two  more  experiments,  in  which  the  rotation  speed  of  the  earth  a >  is  changed 
by  a  factor  of  0.5  (decrease)  and  1.5  (increase).  These  two  experiments  are  referred  to  as 
R45o0.5w  and  R45o  1 .5w,  respectively.  Other  details  of  the  experiments  are  identical  to 
those  of  R45o.  The  oscillation  periods  for  R45o0.5o>  and  R45ol.5a>  arc  16.8  years  and 
37.2  years.  Note  that,  as  <x>  changes,  both  fit  and  /n  change.  Using  the  approximation  of  a 
free  Rossby  wave,  the  phase  velocity  may  be  rewritten  as  c:  cos<p/2/?«  sin:d»,  where  R 
is  the  radius  of  the  earth.  Thus,  as  a>  decreases  (increases),  the  phase  velocity  increases 
(decreases)  and  the  oscillation  period  becomes  shorter  (longer),  consistent  with  the  model 
result.  The  dependence  on  a>  is  not  linear,  because  some  changes  in  c  take  place.  These 
experiments  therefore  lend  support  to  the  Rossby  mechanism. 
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Figure  8.  (a)  and  (b)  Steady  surface  flow  vectors  of  R30o  and  R45o$fatB  respectively.  From  the  steady  state 
of  these  two  runs,  runs  R30p  and  R45ShtBp  are  initiated.  In  each  run  a  heal  anomaly  is  imposed  over  the  three 
northern-most  grid-points  along  the  eastern  boundary,  and  lasts  for  six  months;  (c)  and  (d)  distribution  of  SST 
anomaly  (  C)  5  years  after  the  anomaly  is  removed  for  R30p  and  R45ShtBp  respectively. 


Why  is  the  basin  extent  important  for  the  generation  of  the  oscillation?  In  cases  with 
a  larger  north-south  extent,  maximum  overturning  is  allowed  to  take  place  at  a  higher 
latitude,  and  convection  is  more  concentrated  in  the  north-eastern  comer.  This  favours  the 
generation  of  thermal  w  ind  currents  normal  to  the  weakly  stratified  northern  boundary. 
This  can  be  seen  from  the  steady  surface-flow  of  R30o  and  R45oShtB  shown  in  Figs.  8(a. 
b).  The  flows  of  these  two  runs  are  different.  In  R30o,  in  the  region  near  the  north¬ 
eastern  comer  and  along  the  northern  boundary,  there  are  considerable  northward  flow 
components.  By  contrast,  in  R45oShtB.  there  are  considerable  southward  components. 
The  northward  flows  in  R3l)o  are  important  in  flushing  the  anomaly  in  the  north-eastern 
comer  so  that  it  propagates  along  the  northern  boundary.  (Recall  that  the  adjustment  process 
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in  these  two  runs  is  completely  different.)  This  difference  in  the  flows  accounts  for  the 
different  adjustment  processes  in  the  long-  and  short-basin  cases.  In  R30o.  the  favourable 
flow  condition  leads  to  a  significant  role  of  the  long-Rossby-wave  propagation  at  the  early 
stage  of  the  integration.  This  is  why  strong  oscillations  take  place  (Fig.  3(a)),  although 
they  gradually  weaken  and  are  eventually  damped  by  the  restoring  forcing.  By  contrast,  in 
R45oShtB.  oscillations  arc  absent  from  the  early  stage  of  the  integration  (Fig.  3(b)). 

To  confirm  this  point  further,  two  more  experiments  were  carried  out,  designated 
R30p  and  R45ShtBp  (Table  3)  respectively.  They  started  from  the  steady  state  of  R30o 
and  R45oShtB.  but  each  with  a  heat-flux  anomaly  imposed  for  six  months  in  the  three 
northernmost  grid-points  of  the  eastern  boundary.  At  the  end  of  six  months,  the  heat-flux 
anomaly  was  removed  instantaneously,  and  the  SST  anomaly  in  the  north-eastern  comer 
in  each  case  was  about  5  C.  Figure  8(c.  d)  shows  the  anomaly  distribution  five  years  after 
the  anomaly  was  removed  for  R30p  and  R450ShtBp  respectively.  We  see  that  only  in  R30p 
does  the  anomaly  propagate  along  the  northern  boundary. 

It  should  be  stated  that  the  sensitivity  to  the  model  north-south  basin  extent  may  vary 
with  model  resolution  and  viscosity.  For  reasons  of  numerical  stability,  low-resolution 
OGCMs  usually  require  the  use  of  a  large  viscosity.  It  is  reasonable  to  suspect  that  different 
resolution  and  viscosity  may  lead  to  different  conclusions  with  regard  to  the  effects  of  basin 
extent.  Greatbatch  and  Peterson  ( 1996)  have  shown  that,  in  their  model  under  a  flux  forcing, 
a  smaller  viscosity  leads  to  oscillations  which  are  otherwise  absent.  In  order  to  test  the 
sensitivity  to  viscosity,  we  have  carried  out  tw-o  runs  identical  to  R45oShtB  and  R45o 
except  that  a  horizontal  viscosity  one-third  of  that  in  these  two  runs  was  used.  The  results 
show  thaL  even  with  the  smaller  viscosity,  no  oscillation  is  produced  in  the  short-basin 
case  while  in  the  long-basin  case  persistent  oscillations  occur.  Examination  reveals  that 
the  difference  in  flow  structure  between  the  short-  and  long-basin  cases  discussed  above 
are  seen  in  these  two  new  experiments. 

6.  Discussion  and  Conclusions 

We  describe  an  interdecadal  oscillation  in  a  wind-  and  thermally-driven  OGCM  un¬ 
der  a  restorative  forcing.  Examinations  involving  a  two-dimensional  OGCM,  a  simple 
thermal  ’flip-flop'  model,  and  a  three-dimensional  OGCM  with  or  without  the  nonlinear 
effect  of  temperature  in  the  state  equation  demonstrate  that  the  oscillation  is  not  driven  by- 
mechanisms  apparent  in  the  so-called  convective  oscillator,  or  the  adveciive  overshooting 
oscillator.  Instead,  the  oscillation  occurs  in  association  with  the  generation  of  long  Rossby 
waves  propagating  westwards.  Phase  velocities  predicted  by  theoretical  work  (Gill  1982) 
give  a  lime  for  a  wave  to  travel  from  east  to  west  along  the  northern  boundary  that  is 
approximately  the  same  as  the  oscillation  period. 

That  oceanic  Rossby-waves  play  a  role  in  the  generation  of  decadal  oscillations  has 
already  been  reported  by  other  authors.  Latif  and  Barnett  (1994,  1996)  described  decadal 
modes  over  the  North  Pacific  Ocean  and  North  America  in  a  coupled  ocean-atmosphere 
model.  In  their  model,  oceanic  Rossby-waves  and  their  effects  on  the  large-scale  oceanic 
circulation  largely  determined  the  anomaly  patterns  associated  with  the  decadal  modes. 
The  Rossby  waves  in  their  model  were  driven  by  wind-stress  changes.  The  present  results 
show  that  Rossby  waves  can  be  generated  by  density  changes  and  play  a  significant  role 
in  the  generation  of  decadal  and  interdecadal  oscillations. 

One  of  the  most  interesting  features  of  the  present  study  is  that  the  oscillation  occurs 
under  a  relatively  strong  restoring  boundary  condition.  Winton  (1996)  found  that  a  constant 
heat-flux  forcing  is  required  for  persistent  oscillations.  Under  a  constant  heat-flux  forcing, 
the  zonally  integrated  circulation  of  an  OGCM  is  constrained  by  the  (implied  constant) 
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oceanic  heat  transport.  The  constraint  is  quite  strong,  and  the  imposed  heat-ioss  at  high 
latitudes  must  be  carried  out  even  when  the  model  convection  adjustment  is  switched 
off.  Locally,  the  forced  circulation  is  constrained  bv  a  constant  rate  of  heat  divergence  or 
convergence.  Cai  (1995)  has  discussed  the  consequences  when  these  constraints  cannot  be 
accommodated  by  the  model  internal  circulation.  The  main  finding  in  his  study  is  that  the 
system  will  oscillate  persistently.  The  anomaly  associated  with  the  oscillation  develops 
and  decays  periodically,  there  being  no  damping  from  air-sea  thermal  interactions.  In 
the  present  model,  such  constraints  were  not  applied,  so  that  the  rate  of  convergence  or 
divergence  was  allowed  to  change.  Winton  ( 1996)  reported  a  case  under  a  restoring  forcing 
with  a  restoring  timescale  ten  times  that  used  in  the  spin-up  (equivalent  to  a  restoring  time 
of  250  days  over  a  25  m  mixed  layer),  and  wfith  a  reference  buoyancy -range  increased  by 
a  factor  of  four.  Initial  oscillations  w'ere  seen  but  none  persisted.  In  our  model,  under  a 
30-dav  restoration,  oscillations  were  damped,  but  a  timescale  longer  than  this  produced 
persistent  oscillations.  We  have  shown  that  a  long  north-south  basin  extent  is  conducive  to 
the  generation  of  meridional  flows  normal  to  the  weakly  or  non-strati  fled  boundary'.  These 
flows  arc  crucial  for  the  generation  of  persistent  oscillations.  In  view  of  these  results,  it 
seems  that  the  sensitivity  to  damping  timescale  may  depend  upon  many  model  parameters, 
and  probably  strongly  upon  the  flow  conditions  near  the  polar  boundary. 

The  implication  of  an  oscillation  under  a  constant  heat-flux  forcing  is  that  oceanic 
variabilities  may  be  an  internal  process,  rather  than  forced  externally  as  suggested  by 
Griffies  and  Tziperman  ( 1995).  Such  an  implication  can  be  transferred  to  the  present  study 
because,  firstly,  the  restoring  temperature  is  fixed,  implying  a  fixed  apparent  atmospheric 
temperature  (Haney  1971),  and  secondly,  although  the  heat  flux  is  allowed  to  change, 
the  change  actually  acts  to  weaken  the  oscillation.  However,  in  reality  air-sea  thermal 
interactions  are  present.  The  present  study  shows  that,  even  in  the  presence  of  such  inter¬ 
actions.  oscillations  persist.  This  and  the  role  of  long  Rossby-waves  in  the  generation  of 
interdecadal  oscillations  have  direct  geophysical  relevance. 
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Appendix 

Let  us  simplify  the  "flip-flop*  model  of  Pierce  et  al.  ( 1995)  to  a  thermal  only  system. 
The  salinity  is  set  constant  and  uniform  at  Sn  at  both  layers  and  the  freshwater  flux  is 
set  zero.  The  model  contains  an  upper  and  a  lower  box  with  thicknesses  of  /;.  and  hz, 
and  temperatures  of  7j  and  Tz.  respectively.  Temperature  7)  is  restored  to  T,  at  the  polar 
boundary,  with  a  coupling  strength  of  Kx  equivalent  to  45  days.  As  in  the  OGCM,  surface 
cooling  takes  place.  This  is  associated  w  ith  a  heat  loss,  and,  in  order  to  maintain  a  heat 
budget,  the  lost  heat  is  reintroduced  into  the  system  as  a  subsurface  heating  Q.  The  heating 
can  be  interpreted  as  a  poleward  heat  transport  from  the  subtropics  in  a  more  sophisticated 
model.  In  this  way  Q  is  the  time  average  of  the  heat  loss.  A  vertical  mixing  between  the 
two  boxes  is  allowed  and  the  mixing  coefficient  is  K .  The  coefficient  K  is  taken  to  be 
large  when  the  system  is  statically  unstable,  small  otherwise,  and  is  always  positive. 
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Figure  A.l.  Dependence  of  density  difference  (/>  -  p*)  (kg  m  ')  upon  subsurface  healing  in  a  two-lavcred 
thermal  flip-flop  model.  When  the  subsurface  healing  Q  ranges  from  26  to  45  W  m  the  flip-flop  model  is  in 
an  oscillatory  regime.  In  this  regime,  once  convection  occurs  (large  K )  the  water  column  becomes  stable  (Sf>  <  0). 
whereas  once  convection  stops  (small  K ).  the  water  column  becomes  uastable  {Sp  >  0). 


The  model  equations  for  temperature  are 

hi—  ^K'iT.-TO+Kfc-Tt).  (A.l) 

at 

and 

3  T 

hz—±  =  Q  +  K(Tt  -  Ti).  (A.2) 

dr 

For  simplicity,  the  density  equation  takes  the  form  used  by  Zhang  et  al.  ( 1 992),  i.e. 

p(7\  S)  =  C  +  AS  -  BJ(\  +  B2T).  (A.4) 

In  this  equation,  the  density  (kg  m  3)  is  a  linear  function  of  salinity  5  (parts  per  thousand) 
and  a  nonlinear  function  of  temperature  (  C).  Zhang  et  al.  (1992)  showed  that  when 
C,  A.  £|,  and  Bz  take  values  of  1003.0,  0.77,  0.72.  and  0.72.  respectively,  the  density 
approximates  well  that  in  the  Bryan-Cox  model. 

In  the  steady  state,  the  density  difference  between  the  boxes  is  given  by 

A p  =  p,  -p2  =  K-'<B]Q  +  2B]BzQT,  +  2B]B2Q2/K,  +  B{BZQZ/K).  (A.5) 

A  positive  value  means  that  the  system  is  statically  unstable  and  convection  occurs;  a 
negative  value  indicates  that  the  system  is  statically  stable.  The  system  will  spontaneously 
oscillate  between  convective  and  non-convcctive  states  if  the  steady  state  is  (1)  statically 
stable  (Ap  <  0)  during  convection  (with  a  large  K)  and  (2)  statically  unstable  (Ap  >  0) 
in  the  absence  of  convection  (with  a  small  K). 

Two  features  emerge  from  Eq.  (A.5).  Firstly,  given  a  Q.  the  only  term  that  is  negative 
so  that  Ap  may  be  negative  is  2  B]  BZQT>  (7j  <  0).  This  term  is  the  product  of  the  surface 
cooling,  nonlinearity,  and  the  subsurface  heating.  This  suggests  that  for  oscillation  to 
be  possible,  these  three  processes  must  all  be  present.  In  a  linear  system  Bz  =  0,  Ap 
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is  always  positive,  and  a  non-oscillatorv  convective  state  persists.  We  also  see  that  the 
nonlinearity  acts  to  stabilize  the  water  column,  attempting  to  produce  Ap  <  0.  Further, 
the  only  term  that  can  switch  A p  from  positive  (convective)  to  negative  (non -convective) 
as  K  varies  is  B\B^Ql/K.  This  is  again  dependent  upon  the  nonlinearity.  The  physics 
behind  this  is  simple.  Within  the  water  column,  buoyancy  gain  from  the  subsurface  heating 
provides  the  source  for  buoyancy  loss  by  the  surface  cooling.  If  the  system  is  linear, 
the  buoyancy  gain  from  the  subsurface  heating  is  constant  and  the  system  is  constantly 
convective.  In  the  presence  of  the  nonlinearity,  the  buoyancy  gain  changes  with  T>.  During 
convection,  7)  drops,  and  the  buoyancy  gain  progressively  decreases.  There  is  a  point  at 
which  the  buoyancy  gain  is  too  weak  to  support  convection.  Once  that  occurs  TV  rises  and 
the  buoyancy  gain  progressively  increases.  When  the  destabilizing  buoyancy  gain  is  large 
enough,  convection  is  initiated  again  and  the  cycle  repeats  itself.  Secondly,  for  a  given  /',  if 
Q  is  small  enough,  the  destabilizing  buoyancy  gain  may  not  be  able  to  induce  convection. 
In  this  situation,  the  system  settles  to  a  non-convective  steady  state.  On  the  other  hand, 
if  Q  is  large  enough,  the  destabilizing  buoyancy-flux  is  not  able  to  be  removed  by  the 
given  cooling  strength  and  a  steady  convective  state  results.  These  features  can  be  seen 
in  Fig.  A.l.  which  shows  the  dependence  of  A p  upon  0  given  a  value  of  Ta  that  is  the 
same  as  in  the  polar  boundary  in  R45o  and  a  A'r  equivalent  to  45  days.  Other  values  are 
ht  =  25  m.  B  =  0.072.  A'j  =  6  x  10-7  (when  not  convecting).  and  A?  =  1  x  104  x  A.', 
(during  convection).  W'hen  the  subsurface  heating  Q  ranges  from  26  to  45  W  rrT2.  the 
flip-flop  model  is  in  an  oscillatory  regime. 
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